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ABSTRACT 

In this paper, we take the first step towards an extension of the nonlinear steep- 
est descent method of Deift, Its and Zhou to the case of operator Riemann-Hilbert 
problems. In particular, we provide long range asymptotics for a Fredholm determi- 
nant arising in the computation of the probability of finding a string of n adjacent 
parallel spins up in the antiferromagnetic ground state of the spin 1/2 XXX Heisen- 
berg Chain. Such a determinant can be expressed in terms of the solution of an 
operator Riemann-Hilbert factorization problem. 
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1. INTRODUCTION-MOTIVATION 

It has been established by Korepin, Izergin and Bogohubov [KIB] that many- 
important problems arising in the theory of Quantum Inverse Scattering can be 
reduced to Classical Inverse Scattering Problems. In particular, the correlation 
functions for quantum solvable models can be described by solutions of classical in- 
tegrable equations, or, equivalently, are reducible to the solution of Riemann-Hilbert 
problems. An important observation is that many such quantities are described by 
Fredholm determinants of particular ('integrable') integral operators. As shown 
by different authors (see [DIZ] for a concise proof) the computation of these de- 
terminants reduces to the solution of a Riemann-Hilbert problem. Long range or 
long time asymptotics for such a quantity can then be recovered via an asymptotic 
analysis of an oscillatory Riemann-Hilbert problem.^ In general, such problems are 
not local; in fact they can be considered as operator Riemann-Hilbert problems, 
which makes their analysis harder. In this paper we present a solution of one of 
these problems, by reducing it to two distinct (and essentially different) matrix 
(2x2) problems and thus easily extracting asymptotics. 

We are motivated by the following observation [KIEU]. 

FACT 1.1. Let H(n, t/j, (p) be the probability of finding a string of n adjacent par- 
allel spins up in the antiferromagnetic ground state of the spin 1/2 XXX Heisenberg 
chain. Here ip is an angle related to the external magnetic field /i,^ and 4'{z) is a 
dual quantum field, which can be considered as a holomorphic function in a neigh- 
borhood of the unit circle. Then, n(n, 7/^,0) = ^° det(i+P)'°\ where is the Fock 
vacuum, K is the integral operator on L^[0,oo) with kernel and V is 

an integral operator defined below. 

Our aim is the following. 

THEOREM 1.2. Let ip be as above and (f){z) be an entire function. Let 
P(n, i/j, (j)) — det{I + V), where V is given by (2.1) below. As n — > oo, we have 

(1.1) P{n,'4>,4') ~ const exp(n^log\sin('ip/2) \ + 0(logn)). 

Note that there is no dependence. 

REMARK. Theorem 1.2 falls short of computing the asymptotics for the ac- 
tual correlation function n(n, '0, (/>). This may require the complete asymptotic 
expansion of P{n, </>). We plan to consider this problem in a later publication. 

In the next section, we set up the operator Riemann-Hilbert problem that enables 
us to compute the Fredholm determinant P{n, ip^ (p). Using deformations inspired 
by [DZ2], [FIK] and [K2] we will reduce this problem to two standard matrix ones, 
which we can solve explicitly (as n — > oo). 



^See e.g.[DZl] and [DIZ] for the Ising Chain and the XXO model. Such analysis has also been 
performed in the context of long time asymptotics for soliton equations by several authors; see 
e.g. [DZ2], [Kl], [DVZ], [DKKZ], [K2]. 

^In fact V is defined by e"*'^ = tt < V < 0, where cosh2L = 2/h, L > 0. 
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2. THE OPERATOR RIEMANN-HILBERT PROBLEM 

We are interested in the long range asymptotics of the Fredholm determinant 
det{I + where the integral operator V^") acts on a function / as follows: 

(2.1) (VH/)(^i)= / V^^\z,,Z2)f{z2)dZ2, 

Jc 

where C is the contour 9 ^ z — exp{i9), —ip < 9 < 2tt + ip and 
(2.2) 

V(^){zi,Z2) = -i^ ^[eP(^i)eL")(^2)r(^i,^2)-eL"H^i)eP(^2)r(^2,^i)], 

ZTT Zi — Z2 

Here ip e (— tt, 0) is defined as in Theorem 1.1 and is related to an external magnetic 
field. 

It is established that for special ('integrable') operators where the kernel is of 
the form 

(2.3) V{zi,Z2) = j:^^ii^lM^^ with j:^aj{z)bj{z) = 0, 

zi - Z2 

the computation of the related Fredholm determinant can be reduced to the solution 
of an NxN matrix Riemann-Hilbert problem (see [DIZ] for the proof). In our case, 
the kernel (2.2) is not of the desired type. However, using 

(2.4) r(zi,Z2)= / e 2-1-1 ^ -2--^ ds 

Jo 

we obtain ([FIK]) 

(2.5) V^^\z^,Z2) = -i^ e^f{z,\s)e^::\z2\s)-e^:!\z^\s)e^f{z2\s)ds, 



where 



(2.6) ei") {z\s)^ {z-''expi<Piz) + s^))±V2e-./2. 

z — 1 

In other words, we end up with an operator Riemann-Hilbert problem as follows. 

THEOREM 2.1. [FIK] Let *(z) (similarly M{z)) be a 2x2 matrix of integral 
operators depending on the complex parameter z: 

(2.7) {^jk{z)f){s) = / ^jk{z\s,t)f{t)dy, f e L2[0,oo), 

^0 
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such that 

1. *(2;) is analytic ioi z e C\C. 

2. ^{oo)=l(^ J , where = 5(s -t). 

3. '^+{z) — "if -{z)M{z), z ^ C, where and \E'_ are the normal limits of \E' 
from inside and outside the unit circle respectively, and 

(2.8a) 

Sis -t) + ^e'^^\z\s)e^^\z\t) -7ef )(^|s)eP(^|t) 

je'^^\z\s)e'^I!\z\t) S{s - t) - je^f {z\s)e^::!\z\t) 



M{z\s,t) 
1)1 equivalently, 



(2 8h) M(z) - ( ^ + -Q(^)^--e^(-) , 



where P and Q are given by the kernels 
(2.8c) 



I ^ rl . ^ 2; + 1 S + 

P{z\s, t) = exp[-{s - Y'] 

Q(2;|s,t) = exp[-{s + t)- 



'2' 2 ' 

Then, P(n, V', 0) = (iet(/ + V^"')) satisfies the relation 

Pin +l,ib,d>) 
P(n,7/;,0) 

PROOF: The proof is analogous to the proof of the pure matrix equivalent (see 
[DIZ]); it is essentially given in [FIK]. 

REMARK: We note that P and I — P are orthogonal projection operators. In 
fact, P* =P and P'^{z) = P{z); also P{z)Q{z) = Q{z), and Q{z)Q{l/z) = P{z). 
By definition, P-^(z|s,t) = P{z\t,s). Formulae (2.8c) define Q at least for Rez < 
1. P is defined for at least the unit circle minus the point z = 1. An essential 
singularity exists at 2; = 1. 

From now on, we focus our attention to the (physically interesting case) 7 = 1. 

We use the orthogonal decomposition into Im{I — P) and ImP. Our Riemann- 
Hilbert problem defined by (2.8) splits into two distinct ones, of different nature. 
For convenience, we consider the jump contour as the whole unit circle, with jump 

M* = M, zeC, 
M* = I, ze {\z\ = 1)\C. 

We then have 

(*(/-P))+ = *_(M*(/-P)), 
(*P)+ = *_M*P, 
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across the unit circle, with identity asymptotics at infinity. In other words, we have 
two new problems. Let ^ = ^{I — P) inside the unit circle and ^ = ^ outside the 
unit circle. Then 

(2-10) ,• \ T 

Furthermore, $_|_ has an essential singularity at 2; = 1. 

Likewise, let x — inside the unit circle and x — outside the unit circle. 
Then 

x+ = x-M*P, zeC, 

Again, x+ is essentially singular at 2; = 1. 
Note that 

(2.12) *(0)=x(0) + $(0). 

These new problems are degenerate since none of P, / — P is invertible. They do 
become nondegenerate once we restrict ^{z),x{^) kl < 1 ^^nd M*{z), \z\ = 1 to 
the two orthogonal subspaces I — P and P respectively. Note that for each of the 
new problems the jump across {\z\ = 1}\C is again the identity. Even more, the 
singularity at 2 = 1 is removed. 

On Im{I - P) we have (/ - P)\im{i-p) = I- Thus,^ 

*|/m(/-P),+ = *|/mM(/-P),--^|/m(7-P), where 

(2.13) / -Q(2)2-"e<^(-)\ „ 
M\imii-P)[z) - ^ g(i/^)^ng-0(.) J zee. 

This is essentially the problem appearing in the case < 7 < 1, treated in [FIK]. 
We have the obvious factorization (note Q{z)Q{l/ z)\irn{i-P) = 0) 

M\i^^i_P){z) = Mu{z)Mj^^{z), where 

I -g(2)2-"e^(^) 
/ 



Mu{z) = 

-Q(l/i-e-^W /)• 



By using the analyticity of Q{z) for \z\ < 1 and Q{l/z) for \z\ > 1, we deform the 
problem as follows. 

Let S be the augmented contour consisting of the union of: 



^Note here that Im M{I — P) is independent of z by continuity. 
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1. The contour C. 

2. A smooth curve Cmt joining the endpoints of the contour C, lying entirely 
within the open disc |^| < 1 and close to {|z| = 1} \ C. 

3. A smooth curve Cext joining the endpoints of the contour C, lying entirely in 
\z\ < 1. 

All contours are meant to have a counterclockwise orientation. The complex 
plane is now divided into three regions. 

1. The region containing 0, say Ri. 

2. The unbounded region, say i?2- 

3. The region between Ci^t and Cext, say R^. 

Let 

^ = ^\im{i-P)ML, z e Rl, 

= ^\lTnM{I-P)Mu: Z e R2, 
= ^\lm{I-P) = ^|/mM(7-P), ^ ^ -^3- 

Then, the jumps for $ are equal to Mu on C^xt and on Cint- On C, there 

is no jump since "^^M^^ — ^'_Mi7_. As n — > oo, we trivially get the identity 
solution. So, 

(2.14) X|7m(7-P)(0) ~ /, as oo. 

Note that the essential singularity aX z — 1 never plays a role since z — 1 is not on 
C. 

Also, clearly X|/mp(0) = 0) ^|/m(/-P)(0) = 0- So, it remains to calculate ^|7mp(0); 
in fact, 

(2.15) det*22(0) ~ deW\imP;22{0), as oo. 

Hence, we can simply focus our attention on ImP. 

On /mP, Q{z) = ^^I- In particular, no essential singularity exists and Q{z) 
and Q{l/z) can be defined meromorphically in the whole complex plane using 
Q{z)Q{l/ z) = I. Our operator Riemann-Hilbert problem has jump matrix 

(2.16) M|,„p(.) = /( Q,,/^)^^„^_„., -«Wp*'")..ea 

The conjugation used in the previous case is useless as the resulting scalar problem 
has jump with determinant 0. Recognizing that we have here an operator version 
of a 'shock'-type Riemann-Hilbert problem, we have to use a conjugation with a 
function defined appropriately as a radical (cf. [DVZ], [K2], [DIZ]). 

Let a = —sin'^ip/2 < 0. As in [DIZ], p. 218, we define the function g such that 

g{z) is analytic in C\C; 
g{z) 1, as z — > oo; 

9+{z)g-{z) = -, z e C; 

z 

1^1 <l, z e a 

9- 
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It actually follows that 
(2-17) giz) 



, sog{0) = sin^^. 



We now define a new operator valued 2x2 matrix by 

F^(^) =a^'^^/2$|,^Mp(;2)^(;2)""^a-""^/', \z\ > 1, 



(2.18) 



where as = diag{l, —1) is a Pauli matrix. Then, F„(oo) = / and the jump matrix 
for F„ is 



(2.19) 

As n — > oo, 
(2.20) 



2(fi)"/ -Q{z)e^ 



Q(l/^)e-<^(^) 







zeC. 







1-2" 



1-z 4>(z) 
2 





zeC. 



It is thus appropriate to consider the Riemann-Hilbert 

F(oo) =/. 



(2.21) 



hoping that F = F°° = limn-^ooFn- It turns out that the solution of (2.21) has 
a pole at z = 1. This problem can be conjugated to one with jump equal to 

' i 



-i 



and the solution can be derived easily. We diagonalize M as follows. 



M = SDS-^, 



(2.22) 



S = 



1 



_2i^ -ct>{z) 



2i „-(b(z) J 7 

D = diag{i, —i). 



Now D can be factorized as 

(2.23) 

where 

(2.24) 



D = B_B-^, 
B^diag{p-\P), 



/J(cx>) = 1. 
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Hence {FSB)+ = (FSB)^ across the jump contour; but FSB ~ S near z = oo, 
hence in fact FSB = S and 

(2.25) F-SB S --(^ ^ij3-l3-) 13 + 0-^ 

Even though F has a pole at 2; = 1, the (22)-entry which concerns us does not. It 
is thus possible to show that Fn, 22 — F22, as n 00. The following formulae for 
the Fredholm determinants below follow immediately. 

ip 

detF22{0) = sin—, 

(2.26) ^ 

Using (2.9), (2.17) and (2.18) we get 

%i = .m^"+4(l + 0(^))- 
Pn 2 n 

Theorem 1.2 now follows. 

PROOF OF THEOREM 1.2. A rigorous proof that F°^ = F22 (as weU as 
formula (2.14)) requires a Beals-Coifman type formula which reduces the operator 
Riemann-Hilbert problem to a singular integral equation. This can be easily done 
and the details of the proof are essentially contained in [DIZ], p. 159. The Cauchy 
operator involved takes bounded operators in to bounded operators in and 
is itself bounded. 

A small complication arises because the convergence of is not uniform near 
the endpoints of C. This is dealt with by constructing a parametrix near those 
points, as suggested in [DIZ] for the analogous matrix problem, g being scalar, 
there is no complication due to the operator nature of the underlying problem. 

Finally, one may worry about the fact that F has a pole at z = 1 and the possible 
effect on the validity of the limiting procedure. This, however has been shown not 
to be the case, since any meromorphic problem can be made into a holomorphic 
problem with an extra jump on a circle around the pole (see e.g.[DKKZ]). The 
important fact is that the (22)-entry of F (which is all we are interested in) has no 
singularity. 

The limit F^ is meant to be taken in the trace class norm. The result for 
determinants follows readily and the proof of Theorem 1.2 is complete. 

REMARKS: 

1. It may initially seem that an essential singularity at 2; = 1 is introduced in 
the Riemann-Hilbert problem when one multiples the unknown with the operators 
P and I — P (cf. remark after the proof of theorem 2.1). However, this singularity 
vanishes once we restrict our operators to smaller subspaces. 
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2. In Theorem 1.2, we have been assuming that (f){z) has an analytic extension 
on the complex plane. This is not necessarily true in the original physical problem 
(see Fact 1.1); however, we can always approximate by an analytic that enables 
our calculations to go through, and recover the same result (see [DZ2] for details). 

3. In conclusion, we would like to point out that the importance of the above 
procedure lies in the fact that operator Riemann-Hilbert problems appear in several 
contexts, apart from the computation of correlation functions for exactly solvable 
models of statistical mechanics. For example, it is well known that the inverse prob- 
lem for integrable equations in 2+1 dimensions (Davey-Stewartson, KP, etc.) can 
often be expressed as an operator Riemann-Hilbert problem, and hence deformed 
along the lines described here. 

4. We also note that, even though operators do not generally commute, and 
hence even 'scalar' operator Riemann-Hilbert problems are not trivial to solve ex- 
plicitly, scalar conjugating functions (like g of (2.10)) can still play the important 
'deforming' role they play in the standard matrix factorization problems. 
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